Quantum correlations in a cluster-like system 
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We discuss a cluster-like ID system with triplet interaction. We study the topological properties of 
this system. We find that the degeneracy depends on the topology of the system, and well protected 
against external local perturbations. All these facts show that the system is topologically ordered. 
We also find a string order parameter to characterize the quantum phase transition. Besides, we 
investigate two-site correlations including entanglement, quantum discord and mutual information. 
We study the different divergency behaviour of the correlations. The quantum correlation decays 
exponentially in both topological and magnetic phases, and diverges in reversed power law at the 
critical point. And we find that in TQPT systems, the global difference of topology induced by 
dimension can be reflected in local quantum correlations. 
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I. INTRODUCTION 

Nowadays, quantum correlation has been attracting 
much attention since it plays a crucial role in quantum 
computation and quantum information. Entanglement, 
as an important quantum resource, takes responsibility 
for most quantum information tasks such as quantum 
teleportation and computation But entanglement is 
fragile in open systems. Environment induced decoher- 
ence destroys entanglement correlation in a short time, 
which makes quantum task difficult for implementation. 

However, recent reseach shows that entanglement may 
be not the only worker carrying on quantum tasks. The 
quantum correlation without entanglement may also take 
effect in some scenes, e.g., the quantum computation with 
mixed states plus one pure qubit (DQCl) 0tJ|- 

Quantum discord is developed for the measure of 
"quantumness" of a pairwise correlation [jj. It makes 
it clear that entanglement is one kind of nonclassical cor- 
relation but not all. The quantum discord of some sepa- 
rable states is also nonzero. It may be used as a powerful 
tool to study quantum correlations. 

Lots of work has been devoted to the study of corre- 
lations in different processes, like decoherence and quan- 
tum pha se transition [6l-[l3| . The entanglement of forma- 
tion [14J does not behave smoothly like the correlation 
functions, and shows sudden death and rebirth in some 
scences [15j, which attracts more and more researchers. 
The quantum discord is pointed out to signal the quan- 
tum phase transition Q like fidelity [lj|, while our pre- 
vious work also find that in topological quantum phase 
transition (TQPT) the local correlations are classical and 
the quantum correlation hides in the global system [9j. 

The topological order is a new kind of order beyond 
the conventional symmetry-breaking theory. In topologi- 
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cal order system, the degeneracy of the ground space de- 
pends on the topology of the system configuration, and 
the degenerate ground space is well protected against lo- 
cal perturbation. Such prope rties can be used for fault- 
tolerant computation |17H19| . Another talent idea is 
the measurement-based computation, in which a cluster 
state is prepared and measured as the computation pro- 
cess [20j. There are deep relationship between these two 
methods of computation. 

In this paper, we study the pairwise correlations in a 
ID cluster-like system with triplet interaction, which can 
be implemented in optical lattice [2l|. We discuss the 
properties of the topological order in the system, like the 
boundary dependent degeneracy and topological protec- 
tion. We find the string order parameter (SOP) by the 
method of duality map [22t [23[ to characterize TQPT. 

Furthermore, the system can be decomposed as two 
independent chains of odd and even sites respectively, 
namely, the spin on site i is independent of spin on site 
i + (2n + 1) where n is integer, and we call this "bridge 
correlated". 

The divergency of quantum discord with the distance 
of two site is studied. We find that it behaves in the sim- 
ilar way as the correlation functions, i.e., it decays expo- 
nentially in both topological and magnetic phase areas 
and diverges in reversed power law at the critical points. 

Moreover, the study of the quantum discord and en- 
tanglement shows that the local quantum correlation of 
two sites is suppressed in topological phase area. This is 
different from the study in 2D TQPT [9J, in which local 
quantum correlations vanish completely. And that means 
in TQPT systems the global difference of the topology 
caused by dimension can be reflected in the local quan- 
tum correlations. 

The paper is organized as follows. In Sec. II, we show 
the basic model of this cluster-like system, we discuss the 
topological properties like the degeneracy of the ground 
space and the topological protection. In Sec. Ill, we 
study the quantum correlation of this model. We in- 
vestigate quantm discord, correlation length and mutual 
information. Finally, we draw summary in Sec. V. 
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Figure 1: Configuration of the system in optical lattice imple- 
mentation. Tunneling between the nearest three sites (black 
points) in a triangle gives rise to the triplet interaction term. 



II. TOPOLOGICAL PROPERTIES OF 
CLUSTER-LIKE SYSTEM 

In this section, we introduce the ID cluster-like system 
originally proposed for quantum computation in optical 
lattice [2l| . We calculate the basic property of low-energy 
spectrum. In ID world, there are not too many kinds of 
different topology that we are interested in, except for 
the open string and the closed loop, which correpond to 
open and periodic boundary conditions respectively. We 
show that the degeneracy of the ground space is different 
in these two cases. Besides, the degeneracy is immune 
to local perturbation. These are the typical characters of 
topological order. And we find the SOP to characterize 
the quantum phase transition. 



A. The model of cluster-like system 



whose dimension is just 2 N , i.e., the stabilizers encode 
N — K logical qubits in the protected space. 

Assume we have N sites in the system Eq. (JTJ) , when 
B = 0, it is easy to see that the number of Sj's is N 
under periodic boundary condition, and N — 2 under the 
open boundary condition, the dimensions of the pro- 
tected space of {Si} are 2° = 1 and 2 2 = 4 respectively. 
That is to say, the ground space of the system Eq. ([1]) 
is undegenerated under the loop boundary condtion, and 
4-fold degenerated when it is opened, with the absent of 
external field term. 

Another important property of the system is that all 
the local correlation functions, except those composed of 
products of several Si's, are zero. As an example, for o = 
crfcrj , we can always find certain S k which anti-commutes 
with by so that (6) = (dSk) = — (Sfc6) = — (6) = 0. And 
the system possesses Z 2 -symmetry, i.e., [H, Yl&f] = 0. 
These are important properties as we will see below. 

When B ^ 0, by means of Jordan- Wigner transforma- 
tion 
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we can transform the system Eq. ([T]) into a fermion chain, 
-H/J = ^(c. i _ 1 - Cl t _ 1 )( Cl t +1+Cj+1 ) + S(2cl C4 -l). (3) 



Here, we describe the model we discussed. The Hamil- 
tonian of the system is described as follows, 



Baf 



where J > 0, erf is the Pauli matrix acting on the i-th 
site and Si — <rf_ 1 <T-af +1 . 

The model is originally proposed in Ref. [21] for quan- 
tum computation. It can be implemented in optical lat- 
tice . Atoms are arranged in a triangle lattice as shown 
in Fig. [TJ Tunneling happens in the nearest three sites, 
which gives rise to the triplet interaction term. The one 
body term can be adjusted by Zeeman effect and appro- 
priate laser field. 

When B = 0, the ground space is the common 
eigenspace of Si's that satisfies SJg.s.) = |g.s.), which 
is a typically cluster- like space [241 ]. Cluster states are 
a kind of graph states, which play an essential role in 
measurement-based quantum computation (20| . 

Here we analyze the low ener gy sp ectrum of the system 
Eq. ([TJ) in the stablizer scheme [lj||25|. All Si's commute 
with each other, so we can treat the ground space of the 
system as the protected space of a set of independent 
stabilizer genarators {Si}. 

In the stabilizer scheme, we have N qubits and K in- 
dependent stablizer genarators which are the products of 
Pauli operator erf's. The stablizer genarators commute 
with each other, and the common eigenspace of the sta- 
blizers satisfying Si |$) = | $) is called the protected space. 



We can see that the system can be regarded as two 
independent chains containing odd and even sites respec- 
tively. 

Under the periodic boundary condition, the system 
possesses translational invariance. So it can be diago- 
nalized in Fourier representation, 



-H/J = J2 e^'K - ol fc )(o-fc + 4) + B{2a\a k - 1), 

k 

(4) 

where c„ = ^2 k e lkn ak/VN- By using Bogoliubov trans- 
formation, we get the diagonalized Hamiltonian, 



H/J = 5>(2 7 i 7fc -1), 
k 



(5) 



where e k = (1 + B 2 — 2_Bcos2fe)3, a/. = cos9k~f k + 
ism9kT_ k andtan26>fc = sin 2k/(B — cos 2k). 

When the string is opened, it is difficult to get the low 
energy spectrum and we will discuss the degeneracy of 
the ground space by perturbation method in the following 
part. 



B. Topologically protected degeneracy 

When the string is opened, the ground space is 4-fold 
degenerated when B = 0, as mentioned above. Actually 
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each independent chain contributes two states. In this 
part, we show that this degeneracy is protected against 
external local perturbations. More exactly speaking, the 
energy splitting caused by perturbation tends to zero in 
the thermodynamical limit. 

As the string is opened, Fourier transformation does 
not take effect. We calculate the splitting of the ground 
state energy. Assume the external field is absent at the 
time t — > — oo, and adiabatically switched on. That 
is to say, we construct a new Hamiltonian with time- 
dependent external field A(i) = e' 71 '*, where rj is in- 
finitely small. At t — 0, the system comes back to Eq. 
(P). That is, 



H{t) =H Q + \{t)H'. 



(G) 



Since X(t) is switched on adiabatically, the system evolves 
from the cluster-like ground state |$o) at t — » — oo to an 
eigenstate |G) of Eq. ([1]) at t = 0, which should be one of 
the splitted ground states (2(| [23] ■ The average energy 
of the state |G) is, 



(G|tf(t = 0)|G) 

= ($01^(0, -oo)(H a + H')U(0, -oo)|$ ) ■ (7) 



£7(0,- 



-00 



T exp 



evolution operator, expanded as 



is the time-ordered 



C. String order parameter 

Topological order is an unconventional phase that can- 
not be described by the symmetry-breaking of local or- 
der parameters [28]. When \B\ — y 00, the system leaves 
the topological order and goes to a magnetized phase 
through quantum phase transition. We can find some 
global string order parameter to characterize the phase 
transition. Bellow, we show how to find the SOP by du- 
ality transformation [22|, H3| • 

Under the periodic boundary condition, we make such 
duality transformation bellow to represent the system by 
another self-consistent Pauli algebra {/if}, 



The system turns to be an XF-model, 

-H/J = J2-tftf + i+Btf fJ * +1 . 



(10) 



(11) 



Further, let /j,f 



'i+l 



and 



n 



7-" 

j<i j 



the system 



can be mapped to Ising model in an external field, 



-H/J = J2-t> 



i+l 



'i+2- 



(12) 



,0 

U(0, -00) = 1 + (-i) / dtH a (t) 

J — OO 



(8) 



+ H) 2 / dh dt 2 H' I {t 1 )H' I {t 2 ) + .. 

J — 00 J — 00 

where the perturbation term in interaction picture is 
H' x (t) = \{t)e lHot H'e- lHot 



e -iJt{Si-i+Si +1 ) a z e Ut{Si-x+S i+ i) ^ (g) 



ignoring the boundary terms without loss of generality. 
As e ~ tJtS i = cos Jt — i sin JtSi, we can see that the inner 
product in Eq. (0) is composed of sum of multi-point 
correlation functions, which all vanish until the Ath or- 
der according to what we see in the last part. In the 
Nth term, global terms like (n°f) appear and take ef- 
fect. We can interpret it as a virtual particle running 
along the whole string. Therefore, the energy splitting of 
the ground space is ~ exp(— 1/L), where L is the length 
scale of the system. In thermodynamical limit, L — > 00, 
the degeneracy is perfectly protected like the case in toric 
code [rfj . 

As we see, degeneracy emerges when the loop is 
opened. Besides, the degeneracy is immune against lo- 
cal perturbation when it is not too strong. These prop- 
erties show that the system is a topologically ordered 
system. We can see that there is a close relationship be- 
tween cluster-like system and topological order. Here we 
regard both the topology related degeneracy and topo- 
logical protection as the essential character of topological 
order. 



We can also see that the system is actually composed of 
two independent chains. Combining the two transforma- 
tion together, we can see actually it is 



'i+2) 



'i+l 



'i+l- 



(13) 



The three nearest sites in a triangle (see Fig. [T]) make 
up a new site in the dual lattice. The regular triangles 
and the inverted ones construct two independent Ising 
chains respectively. t£ can be seen as the observable that 
measure the "vortex" of the ith triangle site, clockwised 
or counter-clockwised. 

Lots of work has been devoted to discussing the quan- 
tum phase transition of Ising model. There is a long- 
range order in the dual system (2j|. When \B\ > 1, we 
have 



I™ (t *t* ) = (r£) 2 ~ [1 - 1/B 2 



(14) 



while vanishes when \B\ < 1. rfj can be regarded as 
the order parameter characterizing the phase transition 
at B = ±1. In the original spin representation, we can 
get the hidden SOP as 



A 



even (odd) 



2i(+l)- 



(15) 



Note that we are treating two independent Ising chains. 

When the loop is opened, some boundary terms ap- 
pears, whose effect can be neglected in the thermody- 
namical limit. The physics does not change. 
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Here we emphasize that the existence of SOP is not 
the sufficient condition of topological order. As we see, 
we can also get SOP in the XF-model Eq. (|TT|) . i.e., 

T o T 2j = Yii^Q 1 M? j ^ the duality map, which is a conven- 
sional symmetry-breaking system studied so much. How- 
ever, duality mapping is a useful tool to help us find the 
nonlocal order in topological order system. 



III. PAIRWISE CORRELATIONS 

In this section, we study the pairwise correlations in 
our cluster-like system, like the quantum discord and the 
entanglement of formation (EoF). Quantum discord is 
used as a measure for the "quantumness" of a pairwise 
state. Something interesting are found. We find that 
the quantum correlations are greatly suppressed in the 
topological order area compared with the magnetic po- 
larized phase. The quantum discord decays exponentially 
as the distance of the two spins increases when \B\ ^ 1, 
and diverges in reverse power law at critical points, in the 
behaviour exactly like the two-point correlation function. 
Only the EoF of the spins next-nearest is nontrivial, while 
that of spins further from each other vanishes. 



A. Entanglement, mutual information and 
quantum discord 

Entanglement, as the most important quantum re- 
source, has been discussed lots and there are many dif- 
ferent kinds of measures. One of the most sophisticated 
is the entanglement of formation (EoF) [30|. Entangle- 
ment of formation is an entanglement measure defined 
for bipartite quantum states as 



where H(-) is the Shannon information and H(A\B) is 
the conditional information, which means the average in- 
formation of A we gain when knowing the result of B. 
A natural generalized quantum version is by change the 
Shannon information to von Neumann entropy, 



l(p AB )=S(p A ) + S(p*)-S(p AB ) 



(18) 



Another generalization follows by giving the quantum 
measurement version of conditional entropy. The condi- 
tional entropy implies a measurement on B to get the 
information about A. So we impose projective measure- 
ment {nf } on B and collect the information, 



J(p AB :{IL?}) = S(p A )- 



nfp^nf 



E 



Pi S(il B p AB IL B / Pi ), (19) 



where pi = Tr 

Quantum discord is defined as the minima of the dif- 
ference of X and J [5| , 



DGO=min I{p AB )-J{p AB :{Tlf}) 



(20) 



Due to its power in mixed state quantum computation 
0-51) it has been discussed a lot recently. 

Quantum discord can be used a measurement for the 
"quantumness" of the bipartite correlation. It clears that 
entanglement is not the only "quantum" state. For exam- 
ple, for a separable state p = |00)(00|/2 + | ++)(+ + |/2. 
where |+) = (|0) + \l))/^/2, the quantum discord is not 
zero, which means p contains nonclassical correlation. 



Pairwise state 



E(p) 



mm 



$>s B (|^» 



(16) 



where p is the density matrix of the bipartite states and 
{Pi,\ipi}} satisfies the condition that p — X^J>i|V'i)(V ! i|- 
\tpi) is a bipartite pure state and S (•) gives the von Neu- 
mann entropy of the reduced density matrix of 
For pure states case, this quantity reduces to the entropy 
of entanglement. For two-qubit system, fortunately, EoF 
can be express with concurrence C [l4T |. 

E( P ) = -f(C) log /(C) - (1 - /(C)) log(l - /(C)), (17) 



where /(C) = (1 + Vl - C 2 )/2. The concurrence C = 
max[0, Ai — A2 — A 3 — A4] and A^ are the square roots of 
the eigenvalues of the matrix p(a v (g) a v )p*(a v ® <J y ). 

Mutual information [25] quantifies the amount of com- 
mon information shared by two subsystems. The classical 
mutual information is 

I(A : B) = H(A)+H(B)-H{AB) 
= H(A)-H(A\B), 



To study the correlations in the system, we should first 
get the state of two spins, i.e., their reduced density ma- 
trix. The pairwise density matrix can be decomposed 
by a set of basis {^crf trj}, where and v takes 0, . . . , 3 
and cr? = 1. It can be easily checked that {^cfcrj} 
is orthonormal under the Hilbert-Schmidt inner product 
(A,£)h-s = tr(A^B) H [H[3l|. The reduced density 
matrix of two spins can be written as 



Pij 



4 ^ 



(21) 



where (o>j) = Tr(p G( rfaJ) = trfoyofo?) is the 
Hilbert-Schmidt inner product of p^ and cfcr^. 

Because of the Z 2 -symmetry of the system mentioned 
before, most terms in Eq. (PHI) can be elimited except 
that of lij, cri&j, erf <8> lj and 1, ® er|. So we just 
need to calculate the expectation value of (af c^.^) and 
(<7 Z ) . Since the system can be treated as two independent 
fermion chain like Eq. ([3]), it can be seen that (erf cr^^) 
is zero when R is odd. 
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From the Jordan- Wigner transformation Eq. we 
can get (a z ) and (uq(7p) directly (we take i = without 
loss of generality). 



(cq -cJ)(co+cJ) 
(A B ) , 

((co - 4)(cq + 4)( c fl - C L)( C R + c fi) 

(( c o - cj)(ci + c|)(ci - cj)(c 2 + 4) 
•••(cfl-i -4-i)( c -R + c fi)) 

(-l^-^BoAiBiAa.-.Bji-iAB), 



'o u r/ 



'0 U R/ 



where A, — c,; — cj and £>, = 



. We can check that 



(AoAi) = (BoBj) = when i ^ 0, and the complicated 
expression in the brackets above can be handled with the 
help of Wick theorem [HI Let Gj-i = (AiBj), we 
have 



-R/ 



'0 U BJ 



'0 U R/ 



Go, 

G2 

Go G_i 

Gr-2 Gr-3 

Gi Go 
G 2 Gi 

G/j Gr-\ 



G. 



G-R 

(R-l) 



.. G-i 

G-(R-2) 

G-(r-z) 
G x 



And we have 
G« = 



(A a B R ) = ((c -4)( C fl + 4) 

AT/2 

1 2 e^e^+ > .)(( 7 t-7-p)(7i,+ 



p,q=-N/2 



1 



TV ^ 

P 



e N 



■R„2i8 



N-too 



1 

47T 



2tt 



e 2 



(B - er lr ) 



2tt 



1 + B 2 - 2Bcosr) 1 / 2 ' 



Now we can get the expressions of correlation functions 
above back into Eq. (f/TTj) and we have the reduced density 
matrix of any two spins in the system. Bellow, we will 
discuss the pairwise correlations in the system. 



C. Local correlations in quantum phase transition 

Now we discuss the correlations in the system. First, 
we calculate the EoF of two local spins. As we know, the 
nearest two spins are irrelevant. We give the EoF of the 



o 




Magnet B 

Figure 2: The EoF is shown as a function of the magnetic 
field strength B. The entanglement is born at Be — 0.9767, 
and a sudden change happens at the critical point. 



next-nearest spins shown in Fig. [2] In fact, numerical re- 
sults show that the EoF of the two spins, whose distance 
is further than 2, is zero. 

In Sec. II, we noted that the quantum phase transition 
can be characterized by SOP deduced from duality map- 
ping, and the critical point lies at B = ±1. We can see 
that the EoF in "most" of the topological order area is 
zero and behaves like an order parameter, which is sim- 
ilar to the logarithmic negativity in previous work [2l| . 
However, the EoF is born before reaching B = ±1, at 
the point around \B\ ~ 0.9767. There is a tiny "gap" at 
the critical point, which results from the finite scale, and 
that would eliminate to a singular point in the thermo- 
dynamical limit. 

We cannot treat EoF as an order parameter. However, 
it tells us that in the topological order area, local bipar- 
tite entanglement, as an important quantum correlation, 
is greatly suppressed. It invokes us to study the total 
quantum correlations in this area. 

Second, we calculate the quantum discord of two spins 
with distance R in different magnetic field (Fig. [3]), where 
7q)-8 is even. Around the point B — 1, quantum discord 
has a tiny gap similar to that of EoF. These behaviours 
both root in the property of correlation functions and 
would eliminate to a singular point in the thermodynam- 
ical limit. 

It was mentioned in Ref. [j| that in a 2D TQPT, local 
correlations are always classical and the quantum corre- 
lations hides in the whole lattice, which also happens in 
many other 2D TQPT systems. But things are different 
in ID system (see Fig. [3]). 

In 2D topological order systems, there often exist many 
different conservative string operators whose paths are 
topologically equivalent, and we can always find one that 
anti-commutes with certain local observable. Therefore, 
most local correlation functions would be eliminated and 
the density matrix Eq. (|2ip become diagonalized. 

However, in ID systems, degrees of freedom are re- 
stricted. There are not so many topologically equivalent 
conservative string operators as in 2D. The ID systems 
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Figure 3: The quantum discord is plotted versus B and R where B is the magnetic field strength parameter and R is the site's 
number which correlates with site 0. 




Magnet B Distance 

Figure 4: (Color online) (a) The decay length of quantum correlations (quantum discord, mutual information, and the 
correlation function (<Jq(Jr)) are plotted versus magnetic field strength B. (b) The decay behaviour at B = 1. The correlations 
decay as ~ R~^. We all take R as even. 



do not pocess such high symmetry as in 2D systems, and 
many local correlation functions survive. The quantum 
discord, which measures the quantumness of pairwise cor- 
relations, only gives zero at the cluster state when B = 0. 
Nevertheless, qualitatively speaking, we can see that the 
quantum discord is still quite small in most of the topo- 
logical order area compared with that in the area \ B\ > 1. 
And we say that the local quantum correlation is greatly 
suppressed in the topological phase area. 

On the other hand, this means in TQPT systems the 
global difference of topology induced by dimension is re- 
flected in the local quantum correlations. The dimension 
constrains the topology of the system, and also the types 
of global conservative quantities. In systems with higher 
dimension like Ref. [9], the external field term breaks 
some global conservative operators, while the survival 
ones are still capable to eliminate local quantum correla- 
tions. However, in ID systems like what we study in this 
paper, there are not enough global conservative quanti- 
ties left in th presence the magnetic field and the local 
quantum correlations are just suppressed. The survival 
of local quantum correlation reflects the global restriction 
of the topolgy induced by dimension. 

Besides, we are interested in the decay behaviour of 
quantum discord along with the increase of the distance 
of the two spins we study. Numerical results show that 
the decay behaviours of the quantum discord and total 
mutual information Eq. (|18p are just similar to that of 



two-point correlation functions, i.e., they decay exponen- 
tially when \B\ ^ 1 and with reversed power law at the 
critical points. This is different from the sudden change 
behaviour of EoF, although EoF and quantum discord 
are defined in a similar way, namely, by finding the ex- 
tremum. We show the exponential decay length of the 
correlations with the magnetic field strength B in Fig. [4] 
At the critical points B = ±1, the correlations diverge 
as ~ ij — «, We show them in Fig. HJb). For quantum dis- 
cord, £d ~ 1.0576 and mutual information £m — 1-0179, 
and for the correlation fucntion ((JqOr), £zz — 2.0464. 
We guess this may relate to the universal scaling factor. 
When -6 = 0, the system is the cluster state. All local 
quantum correlations vanish while quantum correlations 
still hide in the chain globally. 

IV. SUMMARY 

We investigate a special model whose Hamiltonian con- 
tains three-spin interactions. This model is composed of 
a cluster and a magnetic term, and we discuss the topo- 
logical properties of this system. The degeneracy of the 
ground space differs in closed and open boundary condi- 
tions, and the degeneracy is topologically protected. We 
obtained the global SOP of this system by the method of 
duality mapping to characterize the TQPT. 

Further, we discuss quantum correlations of this sys- 
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tern. We calculate the quantum discord, mutual infor- 
mation and entanglement in this system. The EoF of 
two local spins is "dead" in most of the topological order 
area. Together with the study of quantum discord, we 
believe the quantum correlation is greatly suppressed in 
the topological order area. This is different from previous 
work in 2D TQPT [9], where local quantum correlations 
all vanish. We believe that is because in ID systems, 
there is not so rich topology or high symmetry as in 2D 
systems. 

On the other hand, in topological order systems, the 
dimension of the configuration constrains the topology of 
global conservative quantities. This global difference of 
topology induced by dimension can be reflected in the 
local quantum correlations. For example, the local quan- 
tum correlations survive in ID TQPT systems, while 
completely vanish in 2D, where there are more global con- 



servative quantities left which root in the richer topology 
of 2D systems. 

Besides, we study the divergency behaviour of the cor- 
relations. The quantum discord and mutual information 
diverge in reversed power law at the critical points and 
exponentially elsewhere. We believe more work can be 
done on the study of the universal scaling behaviour of 
the divergency. 
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